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Abstract 


A  body  of  literature  has  evolved  for  multivariable  systems  which  is  concerned 
with  the  placement  of  closed-loop  eigenvalues  and/or  the  question  of  decoupling.  Atten¬ 
tion  is  turned  to  the  broader  question  of  realizing  specified  rational  transfer  matrices 
with  a  standard  feedback  configuration  for  linear,  time-invariant,  finite  dimensional, 
real,  multivariable,  dynamical  plants  in  this  paper.  A  complete  and  pr  'cise  realiza¬ 
tion  theory  for  asymptotically- stable  plants  is  developed.  Unstable  plants  with  asymp¬ 
totically- stable  hidden  modes  are  also  extensively  treated. 


Introduction 


A  body  of  literature  has  evolved  which  is  concerned  with  the  placement  of  closed- 
loop  eigenvalues  and/or  the  question  of  decoupling  in  multivariable  systems  flj-fll]. 
Attention  is  turned  in  this  paper  to  the  broader  question  of  realizing  specified  rational 
transfer  matrices.  A  standard  feedback  configuration  is  considered  and  attention  is 
restricted  to  linear,  time -invariant,  finite -dimensional,  real,  dynamical  plants.  Spe¬ 
cifically,  the  standard  feedback  configuration  shown  in  Fig.  1  is  studied.  It  is  assumed 


Fig.  1.  Standard  Feedback  Configuration 


that  the  plant,  controller,  and  feedback  network  possess,  respectively,  the  real  state- 
variable  descriptions 


x  =  F  x  +  G  u 

(1) 

~P  P~P  P~P 

y  =  H  x  +  J  u  , 

~P  P~P  P~P 

(2) 

x  =  F  x  +  G  u 
~c  c~c  c~c 

(3) 

y  =  H  x  +  J  u  , 

i.  c  c~c  c~c 

(4) 

and 


h  =  WGf^f 

(5 

(6 

1 


As  a  result  of  the  interconnection  of  plant,  controller,  and  feedback  network  in 
the  configuration  of  Fig.  1  it  follows  that 


Sc  =  S  "  If 


u  =  y 
~p  ic 


2f =  XP  =  l 


The  sizes  of  the  matrices  in  (1)  through  (6)  are  determined  by  the  dimensions  of  the 
vectors  x^,  u^,  y  ,  etc.  The  symbols  used  to  denote  these  dimensions  are: 


v  =  dimx  , 

p  ~p 


V  =  dimx  , 
c  ~c  ’ 


v^.  =  dimx^ 


n  =  dim  y  -  dim  y  =  dim  u,  ,  ( 

m  =  dim  =  dim  yc  ,  ( 

r  =  dim  u  =  dimu  =  dimy,  .  ( 

~  ~c  if  ' 

It  follows  from  (1)  through  (6)  that  the  a  x  m  plant  transfer  matrix  P(e),  the 
m  x  r  controller  transfer  matrix  C(s),  and  the  r  x  n  feedback  network  transfer  ma¬ 
trix  F(s)  are  given  by 

P(s)  =  H  (si  ,  -  F  }* 1  O  +J  ( 

P  vp  P  P  P  ' 

C(s)  =  H  (s  1.  -  F  t"1  G  +  J  ,  ( 

C  v  c  C  C  ' 

c 


F(s)  =  Hf(sl  -  Ff)-1  Gf  +  Jf  .  (18; 

[It  is  evident  from  ( 16 ) -( 1 8)  that  the  symbol  used  in  this  paper  tor  the,  k  x  1  identity 
matrix  is  1^.  ]  The  transfer  matrix  relating  the  transform  of  the  output,  V(s}.  to  the 
transform  of  the  input,  U(s),  is  easily  shown  to  be 


T(s)  =  P(s)Cfs)[ir+F(s)P(s)C(s)]' 


--  [ln+P(s)C(s)F(s)]'J  P(s)C(s)  , 
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Clearly,  (19)  is  meaningful  if  and  only  if  (hereafter  denoted  iff) 

det[lr  +  F(s)P(s)C(s)]  =  det[ln  +  P(s)C(s)F(s)]  *  0  .  (20) 

Moreover,  it  is  desired  that  the  basic  feedback  configuration  be  dynamical.  It  is  shown 
in  [12]  that  this  is  the  case  iff 

lim  T(s)  =  finite  matrix  (21) 

s  -•  « 

or,  equivalently, 

det[lr  +  F(»)P(«)C(“)]  =  det[lr  +  JfJpJc]  i  0  .  (21) 

Practical  arguments  are  also  given  in  [12]  which  justify  the  limitation  of  the  develop¬ 
ments  presented  here  in  accordance  with  the  following  four  restrictions: 

Rj.  The  number  of  plant  inputs  m  equals  or  exceeds  the  number  cf  plant 
outputs  n;  i.  e.  ,  m  >  u. 

R0.  The  number  of  system  inputs  r  equals  the  number  of  system  outputs 
n:  i.  e. ,  r  =  n. 

# 

R^.  The  normal  rank  of  P(s)  is  equal  to  the  number  of  its  rows:  i.e.  , 
normal  rank  P(s)  =  n. 

R^.  The  m  x  n  controller  matrix  C(s)  is  chosen  so  that  the  square  n  x  n 
matrix  P(s)C(s)  has  normal  rank  n:  i.e.,  det[P(s)C(s)]  £  0. 

The  first  significant  contribution  of  the  present  paper  is  best  described  with  the 
aid  of  the  following  definition. 

Definition  1 :  An  n  x  n  rationa'  matrix  T(s)  is  said  to  be  realizable  for  P(s)  if 
for  some  choice  of  asymptotically- stable  dynamical  controller  and  feedback  net¬ 
work  the  standard  feedback  configuration  of  Fig.  1  is  a  dynamical  asymptotically- 
stable  system  possessing  the  transfer  matrix  T(s). 

The  necessary  and  sufficient  conditions  which  T(s)  must  satisfy  in  order  that  it 
v e  realizable  for  P(s)  are  derived  here  for  the  case  in  which  the  plant  is  asymptotically 
stable  and  rank  P(jiu)  =  n  for  all  uj  infinity  included.  It  is  shown  for  this  case  that  the 
limitations  on  the  realizable  T(s)  are  due  to  the  nonminimum  phase  properties  of  the 
plant.  These  properties  are  completely  characterized  for  asymptotically- stable  plants 
by  the  plant  structure  matrix  which  is  introduced  in  the  sequel. 

The  se-ond  significant  contribution  is  the  treatment  of  unstable  plants  whose  un¬ 
controllable  and/or  unobservable  modes  ("hidden  modes",  are  asymptotically  stable. 


*  The  normal  rank  of  a  rational  matrix  is  the  order  of  the  largest  minor  which  is  not 
identically  zero. 
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It  is  shown  that  any  unstable  plant  with  asymptotically- stable  hidden  modes  can  be 
stabilized  with  a  modified  dynamical  observer  of  the  Luenberger  type  f2l.  [13]-[19]. 
Moreover,  the  structure  matrices  of  the  original  and  modified  plant  are  shown  to  be 
strictly  equivalent  and  the  implications  of  this  fact  are  thoroughly  discussed. 

The  notation  used  in  this  paper  is  now  summarized  for  easy  reference,  and  some 
basic  notions  associated  with  a  matrix  function  of  a  complex  variable  "s"  are  defined. 
For  an  arbitrary  matrix  A  the  transpose,  the  complex  conjugate,  the  complex  con¬ 
jugate  transpose,  the  inverse,  the  trace  and  the  determinant  of  A  are  denoted  by 
A  ,  A,  A  ,  A  ,  tr[A],  and  detA,  respectively.  A  diagonal  matrix  A  with  diagona' 
elements  Xj,  . . . ,  X^  is  written  as  A  =  diag[Xj,  X^, ....  X^j.  Column  vectors  are 
represented  by  x,  £,  etc. ,  or  in  the  alternative  fashion  x  =  [x-  x^  •  •  •  xn]f  whenever  it 
is  desirable  to  indicate  the  components  explicitly.  The  n  x  n  identity  matrix,  the  n  x  n 
null  matrix,  the  n-dimensional  zero  vector,  and  the  n  x  m  null  matrix  are  denoted  by 

l,0,o,  and  O  ,  respectively.  The  n-dimensional  column  vector  with  unity 
n  n  ~n  n,  m  (n) 

element  in  the  i'th  row  and  all  other  elements  equal  to  zero  is  denoted  by  e!  or 
simply  e.  when  no  confusion  is  likely  to  result.  The  right  inverse  of  a  p  x  q  matrix 
A  is  the  q  x  p  matrix  A’  which  has  the  property  A  A"  =  1^. 

A  matrix  A(s)  is  rational  when  each  of  its  elements  is  a  rational  function  of  s. 
When  every  element  of  a  rational  matrix  is  finite  at  infinity  it  is  called  a  proper  ma¬ 
trix.  The  matrix  A(s)  is  analytic  in  a  region  of  the  complex  s-plane  when  each  ele¬ 
ment  of  the  matrix  is  analytic  in  the  region.  A  point  s^  is  a  pole  of  A(s)  when  some 
element  of  A(s)  has  a  pole  at  s  =  s^.  A{s)  is  said  to  be  real  if  A(s)  =  A(s).  When  the 
order  of  the  largest  minor  of  A(s)  not  idc  itically  zero  is  v,  then  A(s)  is  said  to  have 
normal  rank  equal  to  V .  Finally,  the  notation 

A^s)  =  A*(-s)  (23) 

is  used  which  for  real  matrices  -  the  only  kind  of  interest  here  -  reduces  to 

A*(s)=A'(-s)  .  (24) 
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Stability  of  the  Standard  Feedback  Configuration 

The  basic  requirements  imposed  on  the  overall  system  in  Fig.  1  are  that  it  be 
dynamical  and  asymptotically  stable.  Conditions  for  the  former  to  be  true  are  stated 
in  the  introduction.  The  latter  requirement  is  discussed  here.  The  first  careful  treat¬ 
ment  of  the  stability  question  for  multivariable  feedback  control  systems  is  due  to 
Chen  [20].  Applications  and  extensions  of  Chen's  results  a  -e  given  by  Youla  [12]. 

Youla  established  for  the  standard  feedback  configuration  of  Fig.  1  the  following  theorem. 

Theorem  1:  When  (22)  is  satisfied,  the  standard  feedback  configuration  is 

asymptotically  stable  iff  the  scalar  function 

Ms)  =  ac(s)Ap(s)Af(s)det[ln+P(s)C(s)F(s)]  (25) 

is  free  of  zeros  in  Res  >  0.  In  (25), 

A  (s)  =  dot(sl  -  F  )  ,  (26) 

c 

Ap(s)  =  det(slv  -  Fp)  ,  (27) 

and 

Af(s)  =  det(slv  -  Ff)  .  (28) 


Iheorem  1  indicates  that  in  general  one  cannot  determine  stability  of  the  standard 
feedback  configuration  solely  from  knowledge  of  the  transfer  matrices  C(s),  P(s),  and 
F(s).  One  must  m  addition  have  knowledge  of  &c(s),  A^fs),  and  A^(s)  which  depend  on 
the  internal  structure  of  the  individual  system  components.  Fortunately,  however, 
practical  considerations  permit  simplifications.  Firstly,  the  controller  and  feedback 
network  are  in  accordance  with  Dex'iait'on  1  to  be  asymptotically  stable.  Thus,  both 
A  (s}  and  A^s)  are  free  of  zeros  in  Rc  s  >  0.  Secondly,  it  is  shown  below  that  one 
can  write 


A  (s)  -  h  (s)ili.  (o) 
P  V  r1 


(29) 


where  hp(s)  is  a  polynomial  whose  zeros  are  associated  with  the  hidden  modes  of  the 
plant  and  (»p{s)  is  the  characteristic  denominator  of  P(s):  i.  e. ,  'ip(s)  is  the  monic 
least  common  multiple  of  the  denominators  of  all  the  minors  of  P(s)  when  these  minors 
are  expressed  as  the  ratio  of  two  relatively  prime  polynomials.  Obviously,  for  every 
practical  plant  h^fs)  is  free  of  zeros  in  Re  s  >  0,  Otherwise,  it  is  not  possible  for  the 
overall  system  to  be  asymptotically  stable.  It  now  immediately  follows  from  Theorem  1 
that 
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Theorem  2;  When  (22)  is  satisfied,  when  the  hidden  modes  of  the  plant  are 
asymptotically  stable,  and  when  the  controller  and  feedback  networks  are 
asymptotically  stable,  then  the  standard  feedback  configuration  is  asymptotically 
stable  iff 

40(s)  =  *  (s)det[l  +P(s)C(s)F(s)]  (30) 

P 

is  free  of  zeros  in  Re  s  >  0. 

Theorem  2  is  significant  in  that  the  test  for  stability  embodied  in  it  can  be 
carried  out  solely  from  knowledge  of  the  transfer  matrices  C(s),  P(s),  and  F(s). 

It  is  now  established  that  (29)  is  a  valid  decomposition.  The  result  follows  from 
the  fact  (see  [21]  and  Theorem  5-19  of  [22])  that  there  exists  a  real  nonsingular  ma- 


where  H  is  n  x  V  ,  G  is  v  x  m,  and 
P  Pi  P  Pi 
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In  addition, 


ys)  =  det(sl^  *Fp) 

and  one  can  easily  establish  from  (31)  that 

A  (s)  =  4{s)detislv  -F  )det(sl 

P2 


or 


h  (s)  =  cet(sl 
P 


v 

P2 


F22)det(slv 


Nonminimum  Phase  Properties 

The  objective  in  this  section  is  the  establishment  of  those  properties  of  the  given 
plant  which  prevent  the  realization  with  F(s)  =  of  an  arbitrarily  specified  rations1 
transfer  matrix  T(s).  These  properties  are  referred  to  as  the  ronminimum  phase 
properties  of  the  plant.  With  F(s)  =  it  immediately  follows  from  Theorem  2  that 
the  overall  system  cannot  bo  asymptotically  stable  vuuess  ^(s)  is  free  of  zeros  in 
Re  s  >  0.  Hence,  closed -right- half- plane  zeros  of  i^fs)  contribute  to  the  nonminimum 
phase  properties  of  the  plant. 

When  F(s)  =  O  then 

T(s)  =  P(s)C(s)  .  (39) 

A  necessary  condition  for  T{s)  -  1r  to  be  realizable  for  P(s)  is,  therefore,  that 

C(s)  =  P_1(s)  ,  (40) 

where  P  *(s)  is  the  ri^ht  inverse  of  P(s).  Since  C(s)  must  be  analytic  in  Re  s  >  0, 
equation  (40)  indicates  that  P  1(s)  must  be  analytic  in  Re  s  >  0.  This  is  never  possible 
when  the  rank  of  F(s)  is  less  than  n,  the  normal  rank  of  P(s),  in  Res  >  0.  Thus,  in 
Re  s  >  0  any  decrease  in  the  rank  of  P(s)  from  its  normal  rank  also  contributes  to  the 
nonminimum  phase  properties  of  the  plant.  The  nonminimum  phase  properties  cited 
ab?ve  are  shown  in  this  paper  to  be  those  properties  of  the  plant  which  restrict  the 
.xass  of  transfer  matrices  that  can  be  realized. 

When  ^(s)  is  free  of  zeros  in  Re  s  >  0,  the  nonminimum  phase  properties  of  the 
plant  are  completely  characterized  by  the  plant  structure  matrix 


r(s) 


That  this  is  the  case  follows  from  the  identity 


■H  (F 
P  P 


"  slv 
_ Pi 

V 

Pi 


•  r(s) 


P(s) 

O 

n,  v 

Pi 

A 

G 

P 

A 

F  -  si 

P  v 

r  n. 

(41) 


(42) 


It  is  clear  from  (36)  that  the  inverse  in  (42)  exists  for  Re  s  >  0  when  £  (s)  is  fret  o* 

P 

zeros  in  Re  s  >  0.  In  this  ca  j,  then,  it  immediately  follows  from  (42)  that 


1 


lhr(s)  =  v. 


pi 


+  rank  P(s),  Re  s  >  0 


(43) 
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Plant  structure  matrices  are  utilized  in  the  sequel  when  the  unstable  plant  is  con¬ 
sidered. 

It  is  possible  to  factor  any  plant  transfer  matrix  into  the  product  of  two  matrices, 
one  of  which  accounts  for  the  nonminimum  phase  properties  of  the  plant.  This  factori¬ 
zation  and  its  properties  are  nor/  discussed.  Given  any  n  x  m  plant  transfer  matrix 
P{s)  of  normal  rank  n  one  ca  i  always  write  provided  the  rank  of  P(ju>)  is  n  for  all 
finite  U)  that 

P(s)  =  V(s)PQ  (s)  „  {44) 

where  the  n  x  m  matrix  Pg(s)  together  with  its  right  inverse  P'^s)  are  analytic  in 
Re  s  >  0  and  the  n  x  n  matrix  V(s)  satisfies 

V*(s)V(s)  =  ln  .  (45) 

The  above  stated  results  follow  easily  when  Theorem  2  of  [21]  is  applied  to  achieve 
the  spectral  factorization 

G(s)  =  P*(s)P(s)  =  PQJ!t(8)P0(s)  .  (46) 

It  is  not  difficult  to  verify  that  the  m  x  m  matrix  G(s)  lias  normal  rank  n  and  that 
the  rank  of  G(jtu)  is  n  for  all  finite  ut).  A  computer  program  for  factoring  G(s)  is 
available  [24]. 

The  paraconjugate  unitary  matrix  V(s)  accounts  for  the  nonminimum  phase 
properties  of  the  plant.  Any  zeros  of  (^(s)  in  Re  s  >  0  are  zeros  of  the  characteristic 
denominator  for  V(s)  and  the  rank  of  V(s)  decreases  in  Re  s  >  0  where  rank  P(s)  does. 
It  is  also  useful  to  note  that  since 

V(s)  =  P(s)Pg1(s)  ,  (47) 

V(s)  is  analytic  in  Re  s  >  a,  o  >  0,  when  P(s)  is  analytic  in  the  same  region.  More¬ 
over,  since  P(s)  is  a  real  matrix  V(s)  is  a  real  matrix,  and  it  follows  from  (45)  that 

V*(juOV(ju>)  =  v'(-jU))V(juu)  =  V*(ju;)V(jU))  =  lu  .  (48) 

It  is  easy  to  infer  from  (48)  that 

lim  V(s)  3  V(®)  =  finite  matrix  .  (49) 

S  —  « 

In  addition  to  the  properties  already  cited  for  V(s)  one  has  from  (47)  and  the  fact  that 
Pq(8)  is  unique  to  within  a  constant  real  orthogonal  matrix  multiplier  on  the  left  that 
V(s)  is  unique  to  within  a  constant  real  orthogonal  matrix  multiplier  on  the  right. 
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(50) 


From  (44),  (45),  and  (49)  it  follows  that 

lim  Pq(s)  =  lim  V"*(s)P(s)  =  lim  VJ{[(s)P(s)  =  finite  matrix  ; 
s-*®  s-*®  s-*® 


The  rational  matrix  P^(s)  is,  therefore,  analytic  at  infinity.  Moreover,  rank  P^(s)  =  n 
for  all  Re  s  >  0.  This  last  property  is  a  consequence  of 


po(8lro1<!>  =  ’«  •  1511 

For  if  rank  Pg(s)  <  n  for  some  s  =  s^,  Re  s^  >  0,  then  (51)  requires  that  P^fs)  have 
a  pole  at  s  =  sQ.  But  this  contradicts  the  analyticity  ot  P~*(s)  in  Re  s  >  0.  The  argu¬ 
ments  just  given  apply  equally  well  at  infinity  provided  the  additional  restriction 
rank  P(jm)  =  n  for  infinite  u>  is  imposed.  In  summar’  ,  Pq(s)  is  analytic  and  rank 
Pg(s)  =  n  in  Re  s  >  0  infinity  included  when  rank  P(juu)  =  n  for  all  uu  infinity  included. 

A  method  for  constructing  a  right  inverse  for  Pq(s)  satisfying 


lim  P  *(s)  =  finite  matr.x 
s  -«®  0 


(52) 


is  now  given.  The  construction  is  accomplished  by  introducing  the  change  of  variable 


s 


1  +  z 
1  -  z 


(53) 


This  transformation  maps  the  region  Res  >  0  of  the  complex  s-plane  into  the  region 
|z|  <  1  in  the  complex  z-plane.  Clearly,  z  =  1  +  jO  is  the  mapping  of  all  points  in  the 
s-plane  infinitely  far  from  the  s-plane  origin.  The  matrix 


W(z) 


(54) 


is  next  considered.  In  view  of  the  properties  of  P^fs),  it  follows  that  W(z)  is  analytic 
in  |z|  <  1  and  rank  W(z)  =  n  for  |z|<  1.  The  matrix  W(z)  therefore  has  the  Smith- 
McMillen  representation  [25],  [26] 

W(z)  =  M(z)[A(z)|On(m_n]N(z)  ,  (55) 

where  M(z)  and  N(z)  are  elementary  polynomial  matrices  of  appropriate  size  and 

A(z)  =  diagfXjfz),  ^2(z) . Xn(z)]  •  (56) 

The  rational  functions  X^(z)  are  all  analytic  in  |z|<  1.  Moreover,  X.(z)  4  0  for  any  z 
satisfying  |z|<  1.  For  if  the  contrary  is  true  then  rank  A(z)  <  n  for  some  z  satisfying 
|z|<  1  which  contradicts  rank  W(z)  =  n  in  the  region  j  z  1.  Hence,  for  any  real 
rational  matrix  K(z)  analytic  in  |z|<  1 
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w'  \z)  =  N-1(z)  [  -^)Z-  ] 


is  a  right  inverse  of  W(z)  which  is  analytic  in  |z|<  1.  It  immediately  follows  that 


P-1!,)  =  »-'(*>  |  ^ 


is  a  right  inverse  of  P^(s)  analytic  in  Re  s  >  0  infinity  included.  The  above  results 
are  summarized  in 


Lemma  1:  A  real  rational  proper  n  x  m  matrix  P(s)  of  normal  rank  n  satisfying 
rank  P(jtu)  =  n  for  all  U)  infinity  included  is  expressible  as  P(s)  =  V(s)Pq(s)  where 
the  n  x  n  matrix  V(s)  and  the  n  x  m  matrix  Pg(s)  are  both  real  rational  proper 
matrices  having  the  properties: 

a)  V*(s)V(s)  =  ln. 

b)  When  P(s)  is  analytic  in  Res  >  a,  a  >  0,  then  V(s)  is  analytic  in  the  same 
region. 

c)  All  zeros  of  the  characteristic  denominator  of  P(s)  in  Re  s  >  0  are  zeros 
of  the  characteristic  denominator  of  V(s). 

d)  The  rank  of  V(s)  decreases  in  Re  s  >  0  wherever  the  rank  of  P(s)  does. 

e)  Both  Pq(s)  and  P~*(s)  are  analytic  in  Re  s  >  0  infinity  included. 

f)  Pq(s)  *s  unique  to  within  a  real  constant  orthogonal  matrix  multiplier  Q 

on  the  left  and  V(s)  is  unique  to  within  the  matrix  multiplier  Q*  on  the  right. 
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The  Main  Theorem  on  Realizability  of  T(s) 


For  the  class  of  asymptotically- stable  dynamical  plants  satisfying  rank  P(ju>)  =  n 
for  all  u)  infinity  included  one  can  always  choose  F(s)  =  and 

C(s)  =  Pg\s)L(s)  ,  (59) 

where  the  n  x  n  real  rational  proper  matrix  L(s)  is  analytic  in  Re  s  >  0  but  is  other¬ 
wise  arbitrary.  For  this  choice  of  controller 

lim  C(s)  =  P0*(m)L(®)  =  finite  matrix  ,  (60) 

s  -*  0 

and  any  minimal  realization  (completely  controllable  and  completely  observable  realiza¬ 
tion)  of  C(s)  is  asymptotically- stable.  Moreover, 

'  (s)  =  P(s)C(s)  =  VfsJPjjWPj^WMs)  =  V(s)L(s)  (61) 

is  realized.  The  traisfer  matrix  T(s)  is  real,  rational,  and  proper;  the  overall  system  - 
plant  with  controller  -  is,  therefore,  dynamical.  Moreover,  V(s)  is  analytic  in 
Re  s  >  0  since  P(s)  is.  Thus,  T(s)  is  analytic  in  Re  s  >  0  and  the  system  is  asymp¬ 
totically  stable. 

The  above  observations  show  that  a  sufficient  condition  for  T(s)  to  be  realizable 
for  P(s)  when  the  plant  is  asymptotically  stable  .a  rank  P(j uu)  =  n  for  all  uu  infinity 
included  is  that  T(s)  =  V(s)L(s),  where  V(s)  and  L(s)  are  as  previously  defined.  It  is 
now  established  that  this  structure  for  T(s)  is  also  necessary. 

Multiplying  (19)  on  the  left  by  F(s)  one  obtains 

F(s)T(s)  =  [lr+F(s)P(s)C(s';- lr][lr  +  F(s)P(s)C(s)r1  (62) 

or 

F(e)T(s)  =  lr.  [lr+F(B)P(s)C(8)]"1  .  ,v3) 

Since  F(s)  and  T(s)  must  both  be  analytic  in  Re  s  >  0,  it  follows  from  (63)  that 
f  1  +  F(s)P(s)C(s)1  1  must  also  be.  Thus, 

T(s)  -  V(s)P0(s)C(s)flr+F(s)P(s)C(s)r1  (64) 

and 

L(s)  r  Po(s)C(s)fJr+F(s)P(s)C(s)r1  (65) 

is  analytic  in  Re  s  >  0.  Moreover,  L(s)  ic  real  and  rational,  and  one  can  establish 
with  the  aid  of  (22)  that  L(s)  is  proper  as  well  when  F(s),  P(s),  and  C(s)  are  real, 
rational,  proper  matrices. 
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The  above  results  are  summarized  in  the  following  theorem: 

Theorem  3:  Given  a  dynamical  plant  with  asymptotically- stable  hidden  modes 

and  a  real,  rational,  proper,  n  x  m  transfer  matrix  P(s)  having  the  properties: 

a)  Normal  rank  of  P(s)  is  n  <  m, 

b)  P(s)  is  analytic  in  Re  s  >  0  infinity  included, 

c)  Rank  P(jUJ)  =  n  for  all  ui  infinity  included,  then  the  necessary  and  sufficient 
condition  for  T(s)  to  be  realizable  for  P(s)  is  that  T(s)  =  V(s)L(s)  where 

d)  L(s)  is  any  real,  rational,  proper,  n  x  n  matrix  analytic  in  Res  >  0, 

e)  V(s)  is  determined  by  the  factorization  P(s)  =  VfsJP^fs)  described  in  Lemma  1. 

Theorem  3  is  the  main  theorem  on  the  realizability  of  T(s).  It  is  restricted  to 
plants  which  satisfy  rank  P(jiu)  =  n  for  all  U)  infinity  included.  For  plants  with  transfer 
matrices  whose  rank  is  less  than  n  at  points  on  the  imaginary  axis  it  is  shown  in  the 
appendix  that  it  is  possible  to  factor  the  plant  transfer  matrix  so  as  to  obtain 


P(s)  =  VTT(s)Pq(s) 


(66) 


where  the  n  x  m  matrix  P  (s)  is  analytic  in  Re  s  >0  infinity  included,  rank  P  (jiu)  =  n 
for  all  U)  infinity  included,  and  Vn(s)  is  analytic  in  Re  s  >  0  infinity  included.  The 
transfer  matrix  P  (s)  can  be  factored  in  accordance  with  Lemma  1  to  obtain 

q 


Pq(s)  =  V(s)P  (s) 


(67) 


where  V(s),  P^fs),  and  P^  (s)  are  analytic  in  Re  s  >  0  infinity  included.  Combining 
(66)  and  (67)  yields 


P(s)  =  VTT(s)V(s)Po(s)  .  (68) 

It  immediately  follows  for  the  choice  F(s)  =  Or  and  C(s)  =  P~*(s)L(s),  where  L(s)  is 
any  real,  rational,  proper  matrix  analytic  in  Re  s  >  0,  that 

T(s)  =  Vn(s)V(s)L(s)  .  (69) 


A  sufficient  condition  for  T(s)  to  be  realisable  for  P(s)  is,  therefore,  that  it  be  factor¬ 
able  in  accordance  with  (69).  This  condition  is  not  necessary,  however.  Other  me¬ 
thods  for  factoring  P(s)  exist  and  the  representation  (68)  is  not  unique. 
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Unstable  Plants 


The  preceding  developments  establish  for  a  large  class  of  asymptotically- 
stable  plants  the  transfer  matrices  T(s)  which  can  be  realized  with  the  standard  feed¬ 
back  configuration.  The  class  of  plants  for  which  P(s)  is  not  analytic  in  Re  s  >  0  is 
treated  in  this  section.  Attention  is  restricted  to  those  plants  with  asymptotically- 
stable  hidden  modes:  only  plants  of  this  type  are  practical.  Before  proceeding,  it  is 
important  to  establish  certain  facts  which  justify  the  procedure  introduced  in  the  se¬ 
quel. 

When  P(s)  is  not  analytic  in  Re  s  >  0,  the  characteristic  denominator  s ) 

contains  zeros  in  Re  s  >  0.  It  immediately  follows  from  Theorem  2,  then,  that  the 
standard  feedback  configuration  cannot  be  asymptotically  stable  with  F(s)  =  O^.  This 
fact  prevents  the  extension  of  Theorem  3  to  unstable  plants. 

A  more  striking  difficulty  with  the  standard  feedback  configuration  is  that  the 
class  of  T(s)  realizable  for  P(s)  can  be  empty  for  some  unstable  plants.  A  simple 
example  is  the  single-input-output  plant  whose  transfer  function  is  (s  -  2)/{s  -  l)(s  -  3). 

It  is  not  difficult  to  establish  that  there  exists  no  asymptotically- stable  controller  and 
feedback  network  which  yields  an  asymptotically- stable  standard  feedback  configuration 
for  this  plant.  This  result  suggests  the  need  for  additional  elements  to  first  stabilize 
the  plant  before  including  it  in  a  standard  feedback  configuration.  In  order  to  handle 
all  cases,  the  additional  elements  should  be  sufficiently  general  so  that  they  permit 
the  stabilization  of  any  unstable  plant  with  asymptotically- stable  hidden  modes.  It  is 
shown  below  that  any  plant  of  the  type  just  described  can  be  stabilized  using  a  modi¬ 
fied  Luenberger  observer  [2],  [13]-[19]. 


In  view  of  (31)  thru  (34),  there  is  no  loss  in  generality  in  assuming  that  the  plant 

has  the  state  variable  description  (1)  and  (2)  in  which  F  ,  G  ,  and  H  are  given,  re- 

I 

spectively,  by  (31)  thru  (33)  with  K  =  lv  .  The  design  of  the  modified  Luenberger 

p 

observer  begins  with  the  formation  of  the  modified  plant  output  vectoi 


where 


y  =  E(y  -  J  u  )  =  EH  x 

ip  'ip  p-p  p~p 


E  =  [e.  e.  •  •  •  e.  1 

L~11~12  ~V 


(70) 


(71) 


,  i^  are  the  numbers  of  the  h  linearly  independent  rows  of  H^.  When 
,  it  follows  from  the  assumption  that  normal  rank  of  P(s)  is  n  that  h  =  n. 


and  ij, 

J  =  O 
p  n,  m 

In  general,  however,  rank  =  h  *  n  is  possible. 

The  plant  state  vector  can  be  written  as 

x  =[x'  x'  x  '  I  , 
~p  [~p1-p2~p3j 


(72) 


1'.- 


where 


dim  x  =  V  ,  i-1,2.3 
-Pi  Pi 

The  objective  is  the  design  of  an  observer  with  state  vector 

z  =  Tx  +  e  , 

-Pi  - 

where  the  error  vector  e  is  exponentially  asymptotically  stable:  i.  e. ,  ||  e  ||  =  Je*e 
<  ce  'v  for  real  constants  c  >  0  and  X  >  0  for  all  initial  error  vectors  at  t  =  tQ. 
The  observer  dimension  is  given  by 


(73) 


(74) 


(75) 


matrix 


That  >  0  is  an  immediate  consequence  of  the  fact  that  the  rank  of  the  n  x  vn 
11^  is  at  most  Vp^.  The  choice  of  the  x  Vp^  matrix  T  is  discussed  in  the  folio  ving 
paragraphs.  Before  proceeding  it  is  first  noted  that 


y  =  [Hjo  | H  1  x  =  H  x  +  H  x 
-p  p  b.v  ~P  p.p,  ,.p3 


(76) 


where 


~  A 

H  =  EH 
P  P 

and 

H  =  EH 
a  a 


(77) 


(78) 


When  T  is  chosen  so  that  the  v_  x  v_  matrix  T  H* 

P1  P1  P 


T*]  is  nonsingular,  then 


(7<>J 


where 


(80) 


is  an  asymptotic  estimate 


lim  x  =  lim  x 
t  -*  00  ~  Pi  t  ->  00  ~  P \ 


of  x  provided  the  error  vector  e,  is  asymptotically  stable. 

rl  ~ 1 

4  «  j  1 _ _  i _ i _ • _ _ •  it.  u _ _  r  _  _  -l 


of  e.  is  deterr  lined  by  the  behavior  with  time  of  e  and  x 
with  K  =  1  immediately  yields 


P3 


The  asymptotic  stability 
Using  (1),  (31),  and  (33) 
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(81) 


x  =  F,,x 
~P3  33-.Pj 

Since  the  hidde'i  modes  of  the  plant  are  asymptotically  stable,  the  eigenvalues  of 
all  have  negative  real  parts.  Tims,  y  is  exponentially  asymptotically  stable. 

The  determination  of  the  behavior  with  time  of  e  requires  more  work.  The 
dynamical  part  of  the  observer  of  interest  is  described  by 


z  =  Az  +  By  +Cu 
~  ~  ~P  ~P 

Substituting  (74)  and  (76)  into  (82)  and  assuming  that  the  matrix  equations 


(82) 


TF  -  AT  =  BH 
P  P 

TG  =  C 
P 


(83) 


are  satisfied  yields 


e  =  A  e  +  (BH  -  TF,  ,)x 
~  -  a  13  ~i 


P3 


(84) 


when  it  is  recognized  that 


x 

~Pj 


=  F  x  + 
P~Pj 


F, ,  x  +  G  u 
13~p3  p~p 


(85) 


Since  Xp^  is  exponentially  asymptotically  stable,  it  follows  from  (84)  that  e  is  ex¬ 
ponentially  asymptotically  stable  whenever  A  has  only  eigenvalues  with  negative  real 

parts.  It  immediately  follows  from  the  fact  that  [F  ,  H  J  is  a  completely-observable 
.  ^  _  P  P 

pair  that  [F^,  H^i  is  also.  Observer  theory  then  guarantees  that  one  cm  always  find 

matrices  A,  B,  C,  and  T  which  satisfy  (83)  and  the  requirement  that  [IT  |T#]  be  non¬ 
singular  and  A  have  only  eigenvalues  with  negative  real  parts. 

It  is  important  to  note  that  the  design  of  the  observer  described  above  depends 

AAA 

only  on  the  matrices  F  ,  G  ,  H  ,  and  J  .  These  matrices  can  be  taken  as  the  ones 

P  P  P  P 

associated  with  any  minimal  realization  of  the  plant  transfer  matrix.  Fortunately, 
algorithms  are  available  for  generating  minimal  realizations  starting  with  the  plant 
transfer  matrix  (see  [22],  Chap.  6).  This  fact  is  important  since  it  shows  that  the 
observer  can  be  designed  from  knowledge  of  only  the  terminal  properties  of  the  plant. 

The  observer  under  consideration  can  be  incorporated  in  a  feedback  loop  around 
the  origir;.!  plant  as  shown  in  Fig.  2.  In  the  figure, 


~  1-1 
H  1 

_E 
T 


(86) 
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and 


This  subsystem  is  referred  to  as  the  modified  plant  in  the  sequel. 


(87) 


Fig.  2.  Modified  Plant 

It  is  now  established  that  one  can  always  chose  the  feedback  matn:.  KQ  so  that 
the  modified  plant  is  asymptotically  stable.  With 

T|  =  [x;  x '  x*  e/]  (88) 

~  ~Pi~p2~P3~ 

it  is  not  difficult  to  verify  that 


"a  a 

F  +C  K. 

P  P  0 

°v  ,  V 

Pi  P2 

F..  +  G  KnL . H 

13  p  0  1  a 

G  KnL, 

p  0  2 

j  = 

F12+GaK0 

F22 

F,,+G  K-L.H 

23  a  0  1  a 

G  K.L, 
a  0  2 

°v  ,  V 

P3  Pj 

°v  ,  V 

P3  P2 

F 

*  33 

V° 

1 - 

o 

<: 

o 

<: 

*— • 

O 

V0’VP2 

BH  -  TF.- 
a  13 

A 

G 

P 


+ 


u 
-  P 


(89) 


(90) 


when  (84)  is  recalled  and  it  is  recognized  that  for  a  properly  designed  observer 

=  Kft(x  +  L ,  H.  x 


u 

~P 


K_x  +  u  =  Kft(x  +  L.H  x  +L_e)  +  u 
0  ~pj  ~p  0  ~pj  l  a  ~pj  Z  ~  ~p 


The  eigenvalues  of  the  coefficient  matrix  of  p  in  (89)  determine  the  stability  of  the 
system.  It  is  not  difficult  to  show  that  the  eigenvalues  of  this  matrix  are  the  eigen- 

A  A 

values  of  A,  ^33’  an<*  *p  +  ^*p^0'  e^8envalues  °f  A  are  the  observer  eigen¬ 

values  which  are  chosen  to  have  negative  real  parts.  The  assumption  that  the  hidden 
modes  of  the  plant  are  asymptotically  stable  is  equivalent  to  all  the  eigenvalues  of  F,, 
and  F-,  having  negative  real  parts.  Since  [F  ,G  }  is  a  completely-controllable  pair, 

P  P  A  A 

it  follows  that  one  can  always  choose  a  Kq  so  that  the  eigenvalues  of  F^+G^Kq  all  have 
negative  real  parts.  An  algorithm  for  choosing  the  matrix  Kq  is  described  in  [27]. 
Hence,  it  is  always  possible  to  make  the  modified  plant  asymptotically  stable.  The 
above  is  summarized  in 


Theorem  4:  Any  real,  linear,  time-invariant,  finite-dimensional,  dynamical 
plant  with  asymptotically- stable  hidden  modes  can  be  stabilized  using  a  suitably 
designed  dynamic  observer. 

Attention  is  now  turned  to  the  computation  of  the  modified-plant  transfer  matrix. 
It  follows  from  (90)  that 


[H  +J  Kj 
P  p  o 


n,  V 


p2 


|H  +J  KftL.H  |j  K.L  Ip  +  T 
'a  p01a’p02~ 


I'~P 


(91) 


Using  the  fact  that  the  inverse  of  a  block  triangular  matrix  with  two  square  blocks  on 
the  diagonal  is  also  a  block  triangular  matrix  of  the  same  form  one  readily  deduces 

A 

from  (89)  and  (91)  that  the  transfer  matrix  relating  Y^(s)  to  U^fs)  is 


P(s)  =  (HptJpK0)(slv 


A  A  1  A 

-  F  -G  K.f  *G  +  J 

p  p  0  p  p 


(92) 


A 

The  matrix  P(s)  is  the  modified-plant  transfer  matrix. 


It  ’mmediately  follows  from  the  identity 


J 

p 

H+JK. 

P  P  o 

J 

P 

_ i 

1 

m 

Ko 

A 

G 

P 

F  +  G  Kft-  si  , 
p  p  0  v 

Pi 

o. 

<o 

_ i 

A 

F  -  si 

P  v 
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r 

O 

■o’ 

*3 

Pi 

that  the  structure  matrices  for  the  original  and  modified  plants  are  strictly  equivalent: 
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ns)  = 


(94) 


J 

p 

H+JL 

p  p  0 

A 

A  A 

G 

F  +  G  Kn-  si 

P 

p  p  0  v 

PlJ 

is  the  structure  matrix  for  the  modified  plant  and 


r(s) 


(95) 


is  the  structure  matrix  for  the  original  plant.  The  modified  plant  is  designed  to  be 
asymptotically  stable;  its  nonminimum  phase  properties,  hence,  are  completely 

A 

characterized  by  the  structure  matrix  r(s).  The  same  is  not  true  for  the  original 
plant  when  P(s)  is  not  analytic  in  Re  s  >  0.  The  set  of  points  in  Re  s  >  0  where  P(s) 
has  no  poles  is  denoted  by  S.  For  all  s  eS  it  readily  follows  from  (42)  and  (93)  th-  u 
(95)  that 


rank P(s)  =  rank T(s)  -  v  =rankT(s)-v  =rankP(s)  .  (96) 

Pi  Pi 

A  A 

Thus,  for  all  s  eS  rank  P(s)  <  n  where  rank  P(s)  <  n.  It  is  also  possible  for  rank  P(s) 
to  be  less  than  n  at  the  points  in  Re  s  >  0  where  P(s)  is  not  analytic.  The  set  of  all 

A 

points  in  Re  s  >  0  associated  with  nonminimum  phase  properties  of  P(s)  is,  therefore 
a  subset  of  the  corresponding  set  for  P(s). 

The  above  results  suggest  that  the  nonmi  'mum  phase  properties  of  the  modified 
plant  are  often  equivalent  to  or  less  severe  than  hose  of  the  original  plant.  One  is 
tempted  to  conclude  from  this  fact  that  the  class  of  T(s)  realizable  for  the  modified 
plant  in  '.he  standard  feedback  configuration  is  equivalent  to  or  larger  than  that  realiza¬ 
ble  with  the  original  plant  even  when  it  is  possible  to  stabilize  the  standard  feedback 
configuration  without  resorting  to  the  use  of  the  modified  Luenberger  observer.  This 
point  has  not  yet  been  rigorously  established,  however. 


Some  additional  observation  ncerning  the  modified  plant  are  now  made.  It  is 
not  difficult  to  verify  using  well  known  matrix  identities  that  (92)  is  equivalent  to 


P(s)  -  (H  +J  Kn)(sl  -F) 

p  p  0  V  p 

pl 


,1a  a  -1a.. 

lG  1  -  Kn(sl  -F  )  *G 

pl  m  0  v  p  pJ 

F  P  pl 


+  J 


(97) 


When  the  bracketed  inverse  in  (97)  is  factored  to  the  right  and  (35)  is  recalled,  one 
easily  obtains  the  relationship 

P(S)  =  P(„)[1  -K0(.l  -F  r'c  r1  .  (98) 

p  Pl 
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Equation  (98)  clearly  places  in  evidence  the  relation„nip  between  P(s)  and  P(s). 
Since  P(s)  is  analytic  in  Re  s  >  0,  the  nonminimum  phase  properties  of  P(s)  are  de¬ 
termined  by  the  points  in  Re  s  >  0  where  all  n-order  minors  of  P(s)  are  zero.  The 
Binet-Cauchy  formula  leads  to  the  fact  that  each  n-order  minor  of  P(s)  is  the  sum 
of  products  of  n-order  minors  of  P(s)  and  [Imp- K0(slv^  -  F  )_1G  ]_1.  It  is  not  easy, 
therefore,  to  relate  in  precise  fashion  the  nonminimum  phase  properties  of  P(s)  and 
P(s)  for  nonsquare  plants. 

Considerable  insight  is  obtainable  for  square  plants.  In  this  case  m  =  n  and  it 
follows  from  (98)  that  1  P 


detP(s)  = 


detP(s) 


det  [1 


K0<slv 


Pi 


A  _  1  A 

F  j  G  ] 
P  P 


(99) 


or 


detP(s) 


detP(s) 


1 


Using  (36)  and 


det<slv  -Fp)  det(slv  -Fp-W 

Pi  Pi  P 


h(s)  =  det(sl  -F  -  G  K-J 
V  p  p  0' 

Pi 


(100) 


(101) 


in  (100)  gives  the  compact  relationship 


a  4’  (s)detP(s) 
P  h(.) 


(102) 


By  design  the  polynomial  h(s)  is  free  of  zeros  in  Re  s  >  0.  Since  P(s)  is  a  rational 
matrix,  it  is  also  true  that 


a  (s) 

detP(s)  =  6^I)  ■ 


(103) 


where  0 _(s)  and  8  (s)  are  polynomials..  Moreover,  B  (s)  divides  the  characteristic 
HP  p 

denominator  i1  (s)  of  P(s): 

P 

,  4’  (») 

Vs)  =  yir  (104) 


is  a  polynomial  in  s.  Substituting  (103)  and  (104)  into  (102)  yields 

A  A-(s)a(s) 

detPfs)  =  — P - E — 

1  h(s) 


(105) 


*  See  page  9  and  equation  (19)  of  reference  [26], 
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The  zeros  of  detP(s)  id  Re  s  >  0  account  for  the  nonminimum  phase  properties 
of  the  modified  plan!.  Since  h(s)  strictly  Hurwitz,  it  follows  that  the  zeros  of 

A 

detP(s)  in  Re  s  >  0  are  the  zero  '  a  (s)  and  a  (s)  in  Re  s  >  0.  The  zeros  of  a  (s)  i*- 

~  p  p  —  p 

As  s  >  0  are  the  points  in  Re  s  >  0  where  rani.  P(s)  <  n.  Any  zeros  of  4  (s)  in  Re  s  >  0 
are  a  result  of  the  fact  that  P(s)  is  not  analytic  in  Re  s  >  0. 

A 

An  example  which  demonstrates  the  generation  of  nonminimum  properties  in  P(s) 
when  P(s)  is  unstable  is  easily  generated.  It  is  not  difficult  to  verify  for 


that 


r  _i_ 

s-1  ■ 

P(S)  = 

s+l 

s+2 

1 

8  +  1 

•  s-1 

s+2  - 

(106) 


d' t' 


a  (s) 


7¥T  =  sV  *  0  ■  0  • 


(107) 


Also, 


Thu  3, 


Ms)  =  (s-  l)(s+l)(s  +  2)  . 

r 


Ms) 

V8)  (s-1|,s+1) 


(108) 


(109) 


and  detP(s)  =  '<  at  s  =  +1. 


A  special  case  of  interest  is  the  single -input- output  plant:  n  =  =1.  In 

this  case, 

d  (s) 

detP(s)  =  P(s)  =-gE 


(110) 


and 


Thus, 


Ms)  =  8  (s) 
P  P 


a  (s) 


detP(s)  =  P(s) 


(111) 


(112) 


and  the  nonminimum  phase  properties  of  the  modified  plant  are  completely  determine 
by  the  zero"  of  the  original  plant  in  Re  s  >  0. 
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Conclusions 


**or  +Iie  class  of  plants  satisfying  the  conditions  of  Theorem  3  the  transfer  ma¬ 
trices  T(s)  realizable  using  the  standard  feedback  configuration  have  been  precisely 
defined.  Much  work  remains  for  this  class  of  plants,  however.  Fundamental  questions 
are  in  need  of  answers.  Given  that  T{s)  is  realizable  for  P(s)  one  can  expect  in  general 
many  combinations  of  C(s)  and  F(s)  which  yield  the  desired  T(s).  Which  of  these  com¬ 
binations  is  best?  One  possibility  is  to  try  to  determine  that  combination  which  mini¬ 
mizes  in  some  sense  the  sensitivity  of  the  system  to  plant  parameter  variations  and/or 
disturbance  inputs.  The  ideas  developed  in  [28]  and  [29]  may  prove  useful  in  this  re¬ 
gard.  Another  possibility  is  to  select  that  controller  and  feedback  network  having  the 
property  that  the  sum  of  the  dimensions  of  the  state  vectors  for  minimal  realizations 
of  both  these  elements  is  a  minimum.  Some  preliminary  results  in  this  regard  are 
contained  in  Chapter  9  of  [22].  Another  possibility  of  course  is  a  compromise  between 
the  two  already  cited. 

With  regard  to  unstable  plants  the  results  reported  here  —  although  extensive  - 
must  be  viewed  as  preliminary  only.  Much  remains  to  be  done.  Suppose  it  is  possible 
tc  stabilize  a  given  P(s)  using  only  the  standard  feedback  configuration.  What  is  the 
class  of  T(s)  realizable  in  this  case  without  resorting  to  the  addition  of  a  modified 
Luenberger  observer?  jl  es  the  addition  of  a  modified  Luenberger  observer  enlarge 
the  class  of  T(s)  realizable  in  this  same  case? 

Finally,  one  can  question  the  sacredness  of  the  scanoird  feedback  configuration. 
This  configuration  is  only  a  special  case  of  the  system  shown  m  Fig.  3.  This  figure 
represents  all  possible  plant  compensation  schemes.  The  connection  network  is 


Fig.  3.  General  Plant  Compensation  Configuration 
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characterized  by  tfte  real,  constant  connection  matrix  M  which  relates  the  intercon¬ 
nection  of  the  syfttem  input  u  and  the  inputs  and  outputs  of  the  compensation  and  plant: 


u 

[fc 

=  M 

y 

u 

£c 

.  u  . 

The  connection  network  can  include  operational  amplifiers  and  this  fact  permits  one  - 
for  all  practical  purposes  -  to  assume  that  the  elements  of  M  can  take  on  any  real 
value.  The  compensation  network  includes  all  dynamical  components  of  the  system 
except  the  plant.  Any  real,  linear,  finite-dimensional,  time-invariant,  asymptotically- 
stable,  dynamic*!  system  .'s  a  possible  choice  for  the  compensation  network.  One  can 
now  raise  all  the  previous  questions  with  regard  to  the  admissable  classes  of  connection 
matrices  and  compensation  networks  just  described  for  the  configuration  shown  in 
Fig.  3.  Work  directed  toward  answering  these  questions  is  presently  under  way. 
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Appendix 


A  method  for  achieving  the  plant  factorization  (66)  is  now  described.  It  is 
assumed  that  P(s)  is  analytic  in  Re  s  >  0  and  rank  P^u^)  <  n.  \  >  Since  P(s)  is  a 
real  matrix,  it  follows  that  rank  P(-jiu^)  <  n. 


The  matrix 


PkM  = 


+ 


A 

s-juik 


+ 


(A.  1) 


is  now  considered  where  the  matrix  A  is  selected  in  accordance  with  the  following 
considerations.  From  the  fact  that  rank  P(jiu^)  <  n,  it  follows  that  there  exists  a  non¬ 
zero  complex  vector  b  such  that 

b^POu^)  =  o^  .  (A.  2) 


The  vector  b  can  be  written  as 

-  =  ~1  + 


(A.  3) 


where  bj  and  b 2  are  real  vectors.  Two  possibilities  exist:  e  ther  the  vectors  b^  and 
b_  are  linearly  dependent  or  they  are  linearly  independent.  The  former  case  is 
considered  first. 


When  the  real  vectors  b.  and  b_  are  linearly  dependent  one  can  write 


b  i  =  Cj  a 


~2  °2~ 


(A.  4) 


where  Cj  and  c^  are  real  scalars  and  the  real  vector  a  satisfies 


Then 


and 


nplies 


The  choice 


||  a  ||  =  Ja.  a  =  ya#  a  =  1 

M  MM  MM 

b  =  (c.  +  jcj  a  =  c  a  ,  c  4-  0 

M  X  4  M  M 

bVP(ja)k)  =  ca#P(ju>k)  =  o^ 

A  =  A  =  a  a* 


(A.  5) 

(A.  6) 

(A.  7) 

(A.  8) 

(A.  9) 
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is  now  considered.  Equation  (A.  1)  becomes 


where 


Clearly, 


Pk(s)  =  V'^sJPts) 


,  [  2s  a  a'  \ 


V"1^)  =  lim  V'^s)  =  lt 
8  —  “ 


(A.  10) 


(A.  11) 


(A.  12) 


Pk(»)  =  lim  Pk(s)  =  P(»)  ; 

8  —  ' 


thus,  both  Pk(s)  and  Vk  (s)  are  proper  matrices.  Using  the  fact  that 


det(l  +AB)  =  det(l  +BA) 
n  m 


(A.  13) 


(A.  13) 


for  an  arbitrary  n  x  m  matrix  A  and  an  arbitrary  m  x  n  matrix  B,  one  easily  obtains 
from  (A.  11)  that 


detVk  (s)  =  det 


(2s  a'a  \  s^  +  2s  +  i 

1+T^T  =  2  T 

,  3  +u,k/  8  +\ 


(A.  14) 


Equation  (A.  14)  establishes  that 


detVk  (s)  =  nonzero  finite  complex  number,  Re  s  >  0,  s  4  ±  .  (A.  15) 


Hence, 


rankP,  (s)  =  rank  P(s),  Re  s  >  0,  s  4  ±juj. 


(A.  16) 


Equation  (A.  16)  is  important  since  it  shows  that  the  set  of  points  in  Re  s  >  0,  s  i  *  jU)k  ■ 
where  rank  Pk(s)  <  n  is  the  same  set  of  points  in  Re  s  >  0,  s  i  ±  jou^f  where  rank  P(s) 


The  next  point  that  one  needs  to  make  is  that  Pk(s)  is  analytic  at  s  =  ±  j This 
fact  follows  immediately  from 


(s  +  ju),  )P,  (s)  =  (s  +  ju),  )V"  (s)P(s)  =0  . 

lS  =  ±jUlk  k  1 S  =  i  j  U)k 


(A.  17) 
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It  is  also  obvious  from  (A.  10)  and  (A.  11)  that  P^(s)  is  analytic  in  Re  s  >  0,  s  ju)^. 
Thus,  Pk(s)  is  analytic  in  Re  s  >  0  infinity  included. 

The  final  property  of  Pk(s)  which  is  investigated  is  rank  Pk(±  This  is  best 

done  by  first  defining  in  lexicographic  order  all  the  corresponding  n'th-order  minors 
of  Pk(s)  and  P(s).  The  j'th  such  minor  of  Pk(s)  and  P(s)  is  denoted  by  A^pj^and  Ap'^, 
respectively.  It  immediately  follows  from  (A.  10),  (A.  14),  and  the  Binet- Cauchy 
formula  (see  [26],  p.  12)  that 


Since  rank  P(±juuk)  <  n,  it  follows  that  for  each  j 

v. 

“k*  “P 


A{p’j)  =  (s2+  «J)  '  A^'j)  .  v.  >  1  , 


(n,  j) 


where  A  L  ’ J ^  0  and  is  finite  at  s  =  ±ju)k.  Thus, 


v.-l 


A(o’j)=  (s2+2s  +  <JU?)(s2+U)5  3  Ain,j) 


(A.  18) 


(A.  19) 


(A.  2?) 


Clearly,  if  for  any  j  it  is  true  that  =  1  then  there  is  one  n'th-order  minor  of  Pk(s) 

which  is  not  zero  at  s  =  ±  juuk  and  rank  Pk(*  ji«k)  =  n.  When  it  is  not  true  that  =  1  for 

any  j,  then  although  rank  P.  (±jou,  )  is  still  less  than  n  one  has  reduced  the  order  of  the 
2  2  ^  “ 

factor  s  +  U/k  in  each  of  the  n'th-order  minors  of  P(s)  by  one.  The  above  process  can 
then  be  repeated  a  finite  number  of  time  -  provided  each  time  that  the  new  vectors  b. 

Ml 

and  b^  are  linearly  dependent  -  until  a  matrix  is  obtained  whose  rank  is  n  at  s  =  ±ju)k< 

Before  considering  the  case  in  which  the  vectors  bj  and  b ^  are  linearly  inde¬ 
pendent,  some  additional  observations  are  now  made.  First,  the  above  developments 
are  easily  applied  to  the  case  rank  P(0)  <  n.  For  this  case,  one  has  immediately  that 
there  exists  a  real  vector  a  satisfying  (A.  8)  with  =  0  since  P(s)  is  a  real  matrix. 

The  final  observation  is  in  regard  to  the  fact  that  (A.  10)  is  not  the  ultimate  relationship 
sought.  One  needs  instead 

P(s)  -  Vk(s)Pk(s)  .  (A.  21) 


It  is  easy  to  verify  that 


Vs> 


2s  a  a' 

_ **  ** _ 

__  -  - 

s  +  2s  +  u), 
k 


(A.  22) 
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satisfies  V^(s)  (s)  =  1R.  Equation  (A.  22)  exposes  the  fact  that  V^(s)  is  a  proper 

matrix  analytic  in  Re  s  >  0. 

Attention  is  now  turned  to  the  case  where  and  b£  are  linearly  independent.  It 
is  first  noted  that  when  the  complex  vector  a  4  o  satisfies 

mt  n 


**P(jV  =  °m  • 


(A.  23) 


satisfies 


b  =  e“-i0a 


b*p(jV  =  o; 


(A.  24) 


(A.  25) 


Moreover,  one  can  choose  a  so  that 


b*b  =  a*a  =  ||aj|Z  =  1 

M  M  M  M  M 

Also,  6  can  always  be  selected  so  that  with  U  a  real  scalar 


(A.  26) 


Thus,  writing 


a'a  =  e’  ^  b#b  =  u  >  0 


a  =  tl  +  j~2 


(A.  27) 


where  a,  and  a,  are  real  vectors,  leads  to 

(»iti"t2i2)+1(&i+*i^  =  u 

Since  u  is  real,  it  immediately  follows  from  (A.  28)  that 

~2~1  +  ~1  ~2  =  Z~l~2  =  ° 


(A.  28) 


(A.  29) 


(A.  30) 


It  is  also  true  that 


*  *  =2i*i  +  *2*2=1 


(A.  31) 


because  of  (A.  29).  Adding  (A.  30)  and  (A.  31)  yields 


~1~1  =  I  f1  + 


(A.  32) 


Subtracting  (A.  30)  from  {A.  31)  gives 


! 

I 


(A.  33) 


Now 


n  n 

U  =  |a7al  =  |  £  a? |<  £  |a.|2  =  |af  =  1  .  (A.  34) 

i  =  1  '  i  =  1 

Hence, 

0  <  u  <  1  .  (A.  35) 


Since  (A.  33)  shows  that  a_  =  o  when  U  =  1,  it  follows  that  a  ’  s  real  in  this  case  and 

~n  •< 

the  results  already  derived  are  applicable.  It  is,  therefore,  assumed  in  the  sequel 
that 


0  <  u  <  1 

The  choice 

.  * 

A  =  aa 

<w  ex 

in  (A.  1)  is  now  considered.  With  this  choice 


where 


Now 


and 


Pk(s)  =  Vk1(s)P(8)  ■ 


(*  -  /  \ 
a  a  aa  \ 

!n  +  s  -  ju)k  +  s  +  jUJk  J 


aa'  =  (aa*) 


Thus,  (A.  39)  becomes 


v;'(si 


1  + 
n 


Zs{ht\  + tzlJ  , 


2  2 

s  +  uu. 
k 


2  ,  2 
s  +  *k 


Clearly, 


=  ii»v;'w  =  i„ 

s  «-•  ® 


(A.  36) 

(A.  37) 

(A.  38) 

(A.  39) 

(A.  40) 

(A. 41) 

(A.  42) 

(A.  43) 
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and 


Pk(-)=  lim  P  (s)  =  P(-)  ; 

s  -*• 


(A.  44) 


r-1 . 


thus,  both  Pk(s)  and  Vk  (s)  are  again  proper  matrices. 

As  before,  the  determinant  of  Vk*(s)  is  of  interest.  Now,  however,  more  work 
is  required  in  order  to  evaluate  this  quantity.  The  computation  is  facilitated  by  con¬ 
structing  the  orthogonal  matrix 


Q  =  f  3i  Sz  "  *  Sj  • 


where 


Si  = 


Si 


and 


Sz  = 


7-|(i+u) 


Sz 


y^d-u) 


(A.  45) 


(A.  46) 


(A.  47) 


That  and  are  orthogonal  is  an  immediate  consequence  of  (A.  29).  It  is  also  clear 
from  (A.  32)  and  (A.  33)  that 


kill  =  llq2ll  = 


(A.  43) 


It  now  follows  that 


Q'V^sJQ  =  1 


n+fi—r)  !  [(^sisi  +(atl)s2Sz]  Car**sn3 

Xs  +UJk/Lq/  J 

(^IdP)313'2"^)  a2,;il [Sl"' 


q  ]  .  (A. 49) 


or 


(A.  50) 
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Since 


detvr^s)  =  dettQ'vr^sJQ]  , 


(A.  51) 


one  easily  bUins  from  (A.  50)  tha' 


s2+2s+tt)2  +  (1  -U2) 


(A.  52) 


detV^(s)- - 

s  +  (i). 


Now  0  <  u  <  1  and  it  follows  from  (A.  52),  therefore,  tnat 

det  V“*(s)  =  nonzero  finite  complex  number.  Re  s  >  0,  a  4  ±  .  (A.  53) 

Hence, 

rank  Pk(s)  =  rank  P(s),  Re  s  >  0,  s  t  ±ji»k  (A.  54) 

just  as  in  the  previous  case  considered. 

From  (A.  38)  and  (A.  42)  it  is  clear  that  Pk(s)  is  analytic  in  Re  s  >  0,  s  i 
It  is  now  established  that  Pk(s)  is  analytic  at  s  =  ±  juu^  as  well.  One  has  from  (A.  42) 
tiiat 


—  .  -1  ,  .  Willi* 

,  )»kIVk  (s)  -  In  +  (si  j*k)  *  (si  j«>kl 


Hence, 


and 


(s-j%)Vk1(s)| 


s  =  +iu). 


<*1L1  +  ~2~2)  +  j(^2tl  -  tllz]  =  It 


(A.  55) 


(A.  56) 


(s  +  jVvk1(s)|  .  +  -i{tzl'i -l)^2]  =  {t?r)  •  (A>57) 

1  a  —  _  "i  111 


's  =  -)Uik 


Since  (A.  23)  also  implies 


avP(juuk)  =  a"P(-juuk)  =  o 
one  readily  concludes  from  (A.  56)  and  (A.  57)  that 


t 

m 


-l. 


(s  +  jwk)Vk  (s)P(s) 


^  0  . 


s  =  ±jwk 


(A.  58) 


(A.  59) 


Thus,  P,  (s)  is  analytic  in  Re  s  >  0  infinity  included. 
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The  final  property  of  P^(«)  which  is  investigated  is  rank  P^(±  jtu^). 
(A.  IS)  one  now  obtains 


+  2s  +  U)£  +  1-u2)^ 

2  1  j 

8  +  u,k  ' 


Instead  of 


(A.  60) 


Since  0  <  u  *  1,  all  of  the  discussion  following  (A.  18)  is  again  applicable.  Moreover, 
now  when  rank  P^(±  ju)^)  <  n  one  is  assured  that  the  process  ca  l  be  repeated. 

The  last  item  requiring  consideration  is  the  computation  of  V^(s).  Clearly, 


Vk(s) 


Q[Q'v: 


,  -1 

(s)Q]  Q' 


(A.  b„ ) 


or  from  (A.  50) 


(A.  62) 


where 

d(s)  sZ  +  2s  +  u)2  +  (1  -  u2)  .  (A.  63) 

Since  Q  is  a  real  finite  orthogonal  matrix  and  since  0  <  u  <  1,  one  has  that  V^s)  is 
analytic  in  Re  s  >  0  infinity  in<  luded. 

One  case  remains  to  be  considered.  It  is  the  one  in  which  rank  P(jw)  <  n.  The 
matrix 

PJs)  =  V^s)  P(s)  (A.  64) 

is  considered  where 

v”*(s)  =  (1  +aa/s)  .  (A.  65) 

®  11  •*  si 

The  real  vector  a  is  one  which  satisfies  ||a||  =  1  and 

M  M 

a'  P(i®)  =  a'P(«)  =  o'  (A.  66) 
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Clearly, 


liin  —  P  (s)  =  lim  aa/P(s)  =  O  .  (A.  67) 

s-- 8  *  ~~  n-m 

Thus,  Pw(s)  is  analytic  at  infinity  or,  equivalently,  P^*)  is  finite.  That  is,  P^s)  is 
a  proper  matrix  analytic  in  Re  s  >  0. 

It  readily  follows  from  (A.  65)  that 

r„t  V”*(s)  =  l  +  s  .  (A.  68) 


Hence,  for  all  finite  s  i  -1  it  follows  that 

rank  P^s)  =  rank  P(s) 


(A.  69) 


Moreover,  in  place  of  (A.  18)  one  now  has 

Ap,j)  =  {s  +  l)A^,j)  .  (A.  70) 

QD 

The  minor  Ap  ’  ^  is  a  rational  function  of  s  and  the  degree  of  the  denominator  poly¬ 
nomial  is  less  than  the  degree  of  the  numerator  polynomial  for  all  j.  Because  of 
(A.  70)  it  is  clear  that  the  difference  between  the  degrees  of  the  denominator  and 
numerator  polynomials  of  A^’ ^  is  one  less  than  the  same  difference  for  i  p' 

Hence,  either  rank  P^®)  =  n  or  the  process  indicated  here  can  be  repeated  a  finite 
number  of  times  until  a  matrix  with  rank  equal  to  n  at  infinity  is  obtained. 


The  final  item  requiring  consideration  is  the  computation  of  Vm(s).  It  is  easy 
to  verify  that 


VJ  8) 


1 

n 


aa's 


s  +  1 


(A.  71) 


Clearly,  V^s)  is  a  proper  matrix  analytic  in  Re  s  >  0. 

Given  a  plant  transfer  matrix  P(s)  there  is  at  most  a  finite  number  of  points  on 
the  imaginary  axis  of  the  complex  s-plane  at  which  rank  P(s)  <  n.  This  is  so  because 
P(s)  is  rational  and  has  normal  rank  n.  Repeated  applications  of  the  factorizations 
described  in  this  appendix  then  leads  to 


P{s) 


(A.  72) 


where  P  (s)  and 
9 


q 

Vs>  =  n  VB) 

k=l  K 


(A.  73) 


are  both  proper  matrices  analytic  in  Re  s  >  0  and  rank  P^(jw)  is  n  for  all  uu  infinity 
included. 
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